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'xf ! ABSTRACT 

m . 

We present multi-frequency observations with the Very Long Baseline Array of the 
\ compact double radio source B2 2050+36. Our observations are at 0.33, 0.61, 1.67, 2.3, 

and 8.4 GHz, with the 0.61 GHz observations forming the third epoch of observation of 

i-C ■ this source at that frequency. At 0.61 GHz, the structure of B2 2050+36 is dominated 

by two components 56 mas apart. Within the uncertainties of the various measurements, 
this separation has remained unchanged for the past 16 years. Any differential image 
wander caused by refractive interstellar scattering is less than 4 mas. Both the lack 
of differential image wander and the frequency dependence of the angular diameter 

. !_h j of B2 2050+36 below 1 GHz indicate that the electron density power spectrum along 

this line of sight has a spectral index near the Kolmogorov value, with a value of 4 
being highly unlikely. We conclude that diffractive scattering dominates along this line 
of sight; results in the literature indicate that this conclusion also holds true for the 
line of sight to the pulsar PSR B2020+28 (8? 7 from B2 2050+36). Comparison of our 
1.67 GHz observations with those obtained 21 years previously place a limit on the 
projected linear separation velocity of the two components of c. 

Subject headings: galaxies: active — galaxies: individual (B2 2050+36) — ISM: struc- 
ture — radio continuum: ISM — scattering — turbulence 



1. Introduction 



Density fluctuations on AU-size scales in the interstellar plasma may produce a host of ob- 
servable effects including refractive intensity scintillations and image wander. The magnitude of 
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these effects depends crucially upon the spectrum of the interstellar density fluctuations. If the 
spectrum is a power law with a spectral index less than 4, these refractive effects will generally 
be masked by diffractive effects from even smaller scale density fluctuations (Romani, Narayan, &; 
Blandford 1986). Conversely, if the spectrum is a power law with a spectral index greater than 4, 
refractive effects will dominate. For many lines of sight the density spectral index appears to be 
close to, but less than, 4 (Armstrong, Rickett, & Spangler 1995), with important exceptions (e.g., 
Hewish, Wolszczan, & Graham 1985; Clegg, Fiedler, & Cordes 1993; Gupta, Rickett, & Lyne 1994; 
Stinebring, Faison, & McKinnon 1996; Rickett, Lyne, & Gupta 1997). 

Complicating the effort to measure refractive effects is their generally long time scale. The 
relevant time scale is given by the time it takes for an AU-sized density fluctuation to drift past 
the observer. For observations of pulsars, whose velocities can exceed 100 km s _1 , the time scale is 
days. For observations of extragalactic sources, the relevant velocity is a combination of the Earth's 
velocity and random interstellar motions and is more likely to be roughly 25 km s^ 1 , meaning that 
the refractive time scale is months to years. 

The source B2 2050+36 is an ideal source for probing refractive effects. It lies on the outskirts of 
the Cygnus superbubble (Bochkarev & Sitnik 1985), a region known to exhibit enhanced interstellar 
scattering (Fey, Spangler, & Mutel 1989; Lazio, Spangler, & Cordes 1990; Fey, Spangler, & Cordes 
1991; Wilkinson, Narayan, & Spencer 1994; Molnar et al. 1995; Spangler k Cordes 1998). VLBI 
observations reveal its milliarcsecond structure to be dominated by two components separated by 
approximately 60 mas, thus a search for refractive image wander requires only relative position 
measurements. 

Fey & Mutel (1993, hereafter FM93) obtained 0.61 GHz VLBI observations of B2 2050+36. 
Combined with previous measurements in the literature, they showed that scattering is moderately 
strong along this line of sight. They were also able to constrain the amount of differential refractive 
image wander between the two components. Their constraint, however, is based on measurements 
of the separation of the components at two epochs separated by four years. If the refractive time 
scale is longer than four years, their constraint is weakened considerably. Assuming that typical 
interstellar velocities along this line of sight are 50-100 km s" 1 , FM93 argued that the refractive 
time scale is comparable to or less than four years. Obviously, if the typical interstellar velocity is 
more like 25 km s _1 , the refractive time scale is much longer. 

We have acquired a third epoch of 0.61 GHz VLBI observations of B2 2050+36, specifically 
to confront the conclusions of FM93 with another epoch of observations. This third epoch is 
12 to 16 yrs after the first two epochs and should allow us to place much better constraints on 
any refractive wander. We have also acquired simultaneous or quasi-simultaneous observations at a 
range of other wavelengths in order to confirm previous measurements of the strength of interstellar 
scattering along this line of sight. 

In §2 we describe our observations and data analysis, in §3 we discuss the implications of our 
observations for both diffractive and refractive interstellar scattering along this line of sight, and in 
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§4 we discuss limits our observations place on the intrinsic structure and kinematics of the source. 
We present our conclusions in §5. 

2. Observations and Data Analysis 

We observed B2 2050+36 with the VLBA for a total of 10 hours on 2000 July 25. Relevant 
observation parameters are summarized in Table 1. 

The observation frequencies were 0.33, 0.61, 1.67, 2.3, and 8.4 GHz; both right and left circular 
polarization were recorded at all frequencies. The 2.3 and 8.4 GHz observations were obtained 
simultaneously in "S/X mode." Simultaneous 0.33 and 0.61 GHz observations also were obtained 
by setting two of the eight baseband converters (BBCs) at each VLBA antenna to record the two 
senses of circular polarization at 0.61 GHz. Because of severe radio frequency interference around 
0.61 GHz, filters restricted the observable bandwidth for each BBC to 4 MHz. The remaining BBCs 
recorded 0.33 GHz. Formatter problems at the Kitt Peak antenna meant that only nine of the ten 
VLBA antennas were correlated. 

The correlated visibilities were processed in the standard fashion within the NRAO Astro- 
nomical Image Processing System (aips). Fringe-fitting intervals ranged from 10 min. at 1.67 GHz 
to 3 min. at 0.33 and 8.4 GHz. At 0.33 and 0.61 GHz, substantial editing of the calibrated visibilities 
was required to remove radio frequency interference. 

In general the calibrated visibilities were then exported to the Caltech difmap program for 
hybrid mapping. To ensure a robust estimate for the separation of the components (§3.2), the 
0.61 GHz observations were imaged using both aips and difmap; the 0.33 GHz observations were 
imaged using only AIPS. 

Figure 1 shows the images of B2 2050+36 at the various frequencies. The observed rms noise 
level in these images ranges from 0.61 mJy beam -1 at 8.42 GHz to 5 mJy beam -1 at 0.61 GHz and 
are summarized in Table 1. The rms noise levels in our images are roughly within a factor of 2-3 
of the expected value. 

Using the images shown in Figure 1 as a guide, we fit the self-calibrated visibility data with 
multiple gaussian component models. The number of gaussian components fitted ranged from 1 
at 0.33 GHz to 7 at 8.42 GHz. These models are summarized in Table 2. Below we shall make 
use of the fitted angular diameters. At the higher frequencies, we take the uncertainty in a fitted 
angular diameter to be one half of the synthesized beam; at the lower frequencies, where a one- or 
two-component fit suffices, we take the uncertainties to be the square roots of the diagonal elements 
of the covariance matrix obtained from a least-squares fit to the visibility data (provided by the 
UVFIT fitting routine within aips); these values are of order 0.5 mas. 
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Fig. 1. — VLBI images of B2 2050+36. (a) 0.33 GHz. The rms noise level is 3.2 mJy beam -1 , and 
the contour levels are 3.2 mJy beam -1 x —2, 3, 5, 7.1, 10, .... The beam is 60.8 mas x 58.5 mas 
at a position angle of 87° and is shown in the lower left, (b) 0.61 GHz. The rms noise level is 
5 mJy beam -1 , and the contour levels are 5 mJy beam -1 x —3, 3, 5, 7.1, 10, .... The beam is 
27 mas x 19 mas at a position angle of 14° and is shown in the lower left, (c) 1.67 GHz. The rms 
noise level is 0.78 mJy beam , and the contour levels are 0.78 mJy beam -1 x -3, 3, 5, 7.1, 10, 

The beam is 7.3 mas x 4.7 mas at a position angle of —7.9° and is shown in the lower left. 

(d) 2.3 GHz. The rms noise level is 1.1 mJy beam -1 , and the contour levels are 1.1 mJy beam -1 
x —3, 3, 5, 7.1, 10, .... The beam is 5.8 mas x 3.6 mas at a position angle of —9.4° and is shown 
in the lower left. 
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Table 1. VLB A Observing Log 







Recorded 


On-source 


Synthesized 


Image 


Frequency 


Bandwidth 


Polarization 


Time 


Beam 


Noise Level 3, 


(GHz) 


(MHz) 




(hr) 


(mas) 


(mJy beam -1 ) 


0.33 


12 


R,L 


2.88 b 


61 x 58 @ 87° 


3.2 


0.61 


4 


R,L 


2.88 b 


27 x 19 @ 14° 


5 


1.67 


32 


R,L 


2.93 


7.3 x 4.7 @ -7.9° 


0.78 


2.27 


16 


R,L 


2.92 c 


5.8 x 3.6 @ -9.4° 


1.1 


8.42 


16 


R,L 


2.92 c 


1.9 x 1.4 @ -10° 


0.61 



a All images are of Stokes I polarization. 

b Observations at 0.33 and 0.61 GHz were conducted simultaneously. 
Observations at 2.27 and 8.42 GHz were conducted simultaneously. 
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Fig. 1.— (Cont.) VLBI images of B2 2050+36. (e) 8.4 GHz. The rms noise level is 
0.61 mJy beam -1 , and the contour levels are 0.61 mJy beam" 1 x —3, 3, 5, 7.1, 10, .... The 
beam is 1.9 mas x 1.4 mas at a position angle of —10° and is shown in the lower left. 
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Table 2. Gaussian Source Models 



V 


S 


r 


8 


a 


b/a 


4> 


Component 


(GHz) 


(Jy) 


(mas) 


( ) 


(mas) 




( ) 




0.33 


2.17 


0.0 




106.2 


0.76 


—35 


E a 


0.61 


2.48 


0.0 




24.1 


0.82 


24 


E 




39 


56 1 


—82 


26.2 


69 


« 


W 

V V 




1 1 


32 6 


« 

o 


47 5 


o no 


67 


E 


1.67 


1.75 


59.3 


—80 


3.5 


0.90 


—49 


w 




0.79 


5.5 


—82 


4.2 


0.74 


-19 


E 




0.99 


1.3 


-132 


4.0 


0.86 


33 


E 




1 06 


3 4 


42 


4 9 


69 


—77 


E 




16 


8 9 


— 106 

1 \J\j 


1 1 8 


32 


—52 


E 


2.27 


1.97 


0.0 


—32 


2.9 


0.71 


—61 


W 




0.92 


58.5 


101 


3.7 


0.81 


17 


E 




0.58 


53.8 


100 


3.3 


0.59 


-25 


E 




0.83 


61.5 


97 


4.2 


0.54 


-78 


E 




0.12 


52.9 


103 


7.5 


0.39 


-46 


E 


8.42 


0.35 


18.4 


-89 


0.8 


0.54 


84 


W 




0.39 


20.1 


-87 


1.0 


0.51 


-1 


W 




0.19 


34.9 


103 


2.1 


0.64 


-29 


E 




0.19 


40.3 


105 


2.1 


0.54 


60 


E 




0.19 


43.2 


99 


3.5 


0.51 


-85 


E 




0.22 


18.8 


-85 


3.1 


0.24 


-67 


W 




0.16 


39.0 


102.4 


6.3 


0.28 


40 


E 



Note. — The models fitted to the visibility data are Gaussians with 
flux density S and FWHM major axis a and minor axis b, with major 
axis in position angle <p (measured north through east). Components 
are separated from the (arbitrary) origin of the image by a distance r in 
position angle 0, which is the position angle (measured north through 
east) of a line joining the components with the origin. Components 
and sub-components are identified as with either an "E" or "W" corre- 
sponding to whether they are Eastern or Western, respectively. 

a Assumed Eastern component; see §3.1. 
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3. Constraints on Interstellar Scattering 

We now revisit the scattering properties of the line of sight to B2 2050+36. We consider both 
diffractive effects, namely angular broadening, and refractive effects, namely differential image 
wander and refractive intensity fluctuations. 

We shall make the standard assumption throughout that the density fluctuations responsible 
for interstellar scattering can be parameterized by an exponentially-truncated power law 

P Sne (z,q)=C 2 n (z) q - a e-^ 2 , q>q . (1) 

Here q± is the largest wavenumber on which density fluctuations occur, which is related to the 
smallest length scales for density fluctuations l\ or "inner scale" as qi = 2ir/l\. The smallest 
wavenumber on which density fluctuations occur or the "outer scale" is go = 2tt/Iq. The con- 
stant C\ is a measure of the level of density fluctuations and is related to the rms density as 
(nl) = [2(2vr) 4 " a /(a - 3)] C£ ll~ a (Cordes et al. 1991). The power law spectral index is assumed 
to provide clues about the process(es) that generate or maintain the density fluctuations. Our 
observations allow us to measure a, constrain C%, and indirectly place limits on l\. 



3.1. Diffractive Scattering 

3.1.1. Angular Broadening and Density Spectrum 

An infinitely-distant point source viewed through a medium filled with density fluctuations 
having a spectrum described by equation (1) and a spectral index a = 11/3 will be broadened to 
an angular diameter (Cordes & Lazio 1991; Taylor & Cordes 1993) 



f 128 mas SM 3 / 5 ^ 1 / 5 6 d < CTOSS ; 

132 mas SM 1 /^ (_h_\ ~ 6d > 



e d 

Here fGHz is the observing frequency measured in GHz, and the scattering measure SM is the line 
of sight integral of C^, 

SM= / dzCl(z). (3) 



The quantity # C ross is the cross-over angle and describes the relative importance of the inner scale l±. 
It is given by Across — f/ / 16^Qg z (Zi/100km) _1 . The constants (128 mas and 132 mas) given in 
equation (2) are those appropriate for infinitely distant sources and differ slightly from those given 
by Cordes & Lazio (1991) who considered sources embedded in the medium. These formulae specify 
the diameter as the full width at half maximum (FWHM) of the shape of the image and that shall 
be the convention we adopt in reporting the measured diameters. 

The scattering diameter 9d defines the diffractive scale Id, Id = A/#d (at observing wave- 
length A). If the typical interferometer baseline b is sufficient to resolve 6 d (i.e., b ~ l d ) and b <h, 
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we expect Od > 9, 



and Od oc v 2 . Conversely if 6 > h, then we expect 9d < 9, 



and 0^ oc v 



-11/5 



cross 



cross 



Therefore fitting the observed diameter as a function of frequency can constrain both l\ and a. 

There are alternate causes for a 9& oc dependence. For instance, if the spectrum of density 
fluctuations is not a power law but a gaussian centered on a particular length scale, that can produce 
a 9d oc v dependence. Similarly, if the density spectrum is dominated just slightly by AU-size 
scales such that a > 4, that could also produce an apparent Od oc v~ 2 dependence within the 
measurement uncertainties. As we shall show below, a < 4 so that we could interpret a Od oc v~ 2 
dependence as resulting from a "large" inner scale. 

Figure 2 shows the measured angular diameters as a function of observing frequency. Also 
shown are two curves, one for oc v~ 2 - 2 appropriate for a spectrum with a = 11/3 and a diffractive 
scale Id much larger than l\ and one for oc v~ 2 appropriate for the case in which Id <h- 

As Figure 1 demonstrates, B2 2050+36 is not only not a point source, but the two components 
of it can be resolved into multiple sub-components at high frequencies. A visual inspection of 
Figure 2 suggests that at frequencies below 1 GHz scattering dominates the structure of B2 2050+36 
while at frequencies above 1 GHz there are clear deviations from either curve due presumably to 
intrinsic structure. In order to quantify this result, we have fit the measured diameters of the 
eastern component (Table 2) with the following functional form 



Here 0i is the intrinsic diameter and S is the scattering diameter, both at the fiducial frequency 
of 1 GHz. Unlike in equation (2) where we specified the dependence for the specific case of a = 11/3, 
here we allow the frequency dependence to be part of the fit. The frequency dependence we use is 
appropriate for a < 4. 

Several points are in order regarding the choice of data we used. At 0.33 GHz only one 
component is seen (Figure la). Attempts to fit the visibility data with a two-component structure 
suggested by the 0.61 GHz image produced a poorer fit to the data. We have subtracted a one- 
component model (Table 2) from the visibility data and imaged the residuals. There is no indication 
of a second component at the location expected from the higher frequency images. We have also 
produced model visibility data, comprised of two components whose flux densities, diameters, and 
positions are consistent with those seen at higher frequencies. Fitting these model data with a one- 
component model results in a fit nearly the same as what we find by fitting the observed visibility 
data. The similarity suggests that any systematic error in determining the diameter is probably 
no more than 1 mas. Because of the self-calibration we performed, we have no absolute position 
information. Given that the eastern component is stronger at 0.61 GHz, we have therefore assumed 
that the single component seen at 0.33 GHz is the eastern component. 

At the higher frequencies both components can be resolved into multiple sub-components. In 
Figure 2 we plot the diameters of the individual sub-components. The diameters of the individual 
sub-components are clearly not representative of the diameter of the entire component. Therefore, 



2 = 2 + 2 s v- : 




(4) 
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Frequency (GHz) 




Fig. 2. — The angular diameter of the components of B2 2050+36 as a function of observing 
frequency. All components listed in Table 2 are plotted. Top: The eastern component. Bottom: 
The western component. For the western component we also show, as an open circle, the diameter 
it would have assuming that its scattering diameter is identical to that of the eastern component. 
This point is shown for display purposes only and is not used in any of the fits described. In both 
panels a dashed line shows the dependence 6 oc z/~ 2 - 2 , as expected for a Kolmogorov spectrum, 
and a dotted line shows 6 oc v~ 2 . Both lines are normalized to pass through the measured angular 
diameter at 0.61 GHz. For the eastern component only, a solid line shows the fit, which allows 
for intrinsic structure and which is discussed in the text. At the higher frequencies, we plot all 
sub-components, but we show the uncertainty on only one of them as we assume the uncertainty 
to be the same for all sub-components and equal to one half of the synthesized beam. 
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as a measure of the size of the components at these frequencies, we have added the diameters of 
the sub-components (Table 2) in quadrature. If anything, this estimate will underestimate the size 
of a component. However, our primary interest is not in the intrinsic structure, and this procedure 
suffices to obtain an estimate for the intrinsic size. 

B2 2050+36 is classified as a gigahertz-peaked spectrum source, and both components show 
a single flux density maximum near 1 GHz (Mutel, Hodges, & Phillips 1985). Accordingly, in 
equation (4) we have assumed that the intrinsic structure has no frequency dependence. While we 
have not conducted extensive tests, assuming no frequency dependence produces clearly much better 
fits than does assuming a frequency dependence of v~ x , as might be expected for a homogeneous 
source with a peak brightness temperature Tg (Kellermann & Owen 1988). 

We have used a grid-search technique to minimize the difference, in a x 2 sense, between the data 
and the expression of equation (4) for the three parameters 0j, S , and a. We have also determined 
the allowed ranges for the parameters in the following manner: Holding two parameters fixed at 
the values that minimize x 2 ■, we change the third parameter until x 2 changes by unity. 

Table 3 summarizes the best fit value for the parameters and their allowed ranges. The best- 
fitting curve is also shown in Figure 2. A key result of this fit is that a density spectrum spectral 
index of 4 is highly unlikely. This result means that the inner scale l\ is substantially smaller than 
most baselines in the array or that 8d < cross . The typical array baseline is roughly 2000 km, so 
we place an upper limit of l\ <C 2000 km. This constraint is not a serious challenge to the value 
estimated by Spangler & Gwinn (1988) of l\ < 200 km. Using equation (2) we can also derive the 
scattering measure SM. Table 3 also lists the values of these derived quantities. 

3.1.2. Image Anisotropy and Orientation 

The anisotropy and orientation of a scattered image can, in principle, yield information about 
the density structures responsible for the angular broadening (e.g., Wilkinson et al. 1994; Trotter, 
Moran, & Rodriguez 1998; Cordes & Lazio 2001). The simultaneous low frequency observations we 
report here are seemingly ideal for constraining any wavelength-dependent changes in the anisotropy 
or orientation of B2 2050+36, as the u-v plane coverage should be identical (once scaled by the 
ratio of the wavelengths) . 

We consider first the anisotropy of the image. We find axial ratios of b/a ~ 0.80 (Table 2) 
or, in terms of the image anisotropy parameter of Romani et al. (1986), e s ~ 0.25. Romani et 
al. (1986) also provide an estimate of the image anisotropy to be expected from random fluctua- 
tions within the medium. For our observing wavelengths and scattering parameters derived above, 
the expected random contribution to e s is approximately 0.03. We conclude that these image 
anisotropies are significant. These anisotropics are also comparable to those of other sources seen 
through the Cygnus region (Spangler & Cordes 1988; Wilkinson et al. 1994; Molnar et al. 1995; 
Spangler & Cordes 1998; Desai & Fey 2001). Desai & Fey (2001) attribute the similarity of the 
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axial ratios, between different sources at different frequencies and different epochs, to be indicative 
of an anisotropy in the scattering material itself. 

We cannot conclude that the image anisotropy changes with frequency for B2 2050+36. The 
formal uncertainties from the fitting of the visibility data are approximately 1%. Even this small 
uncertainty on the angular diameters, when propagated to the axial ratios, implies an uncertainty 
of A(b/a) ~ 0.03, comparable to the difference between the axial ratios at the two frequencies 
(Table 2). 

Changes in the image orientation between the two frequencies may reflect changes in the 
orientation of the density irregularities responsible for the scattering. Unfortunately, we also cannot 
draw any conclusions about the image orientation. Although the u-v plane coverage should be 
similar, differences in the editing (due to the significant levels of RFI encountered) change the 
relative u-v plane coverage significantly. An indication of the difference can be seen by inspecting 
the synthesized beams shown in Figures la and lb. In particular, the synthesized beam at 610 MHz 
has an orientation of 14° compared to an orientation for the eastern component of 24°. Thus, the 
apparent change in orientation between 330 MHz (—35°) and 610 MHz (24°) may be due largely 
to the differences in u-v coverage. 

3.2. Refractive Interstellar Scattering 

Large scale density fluctuations may affect the image position and flux density. In this section 
we place constraints on the strength of refractive modulation of the image, focussing primarily on 
the 0.61 MHz image. 

3.2.1. Differential Image Wander 

Density fluctuations on a scale l r = DO^, where D is the characteristic distance to the scattering 
medium, will cause refractive "bending" and "steering" of the propagating rays. In turn, this will 
affect the angle of arrival of rays and make radio sources wander about a true or nominal position. 
The magnitude of this effect depends both upon the strength of scattering and the value of a. 

FM93 selected B2 2050+36 because its two component structure allowed for relative position 
measurements to be performed, obviating the need for determining absolute positions with milliarc- 
second accuracy. FM93 summarized the predictions for the amount of differential angular wander 
as a function of a as predicted by various groups (Cordes, Pidwerbetsky, & Lovelace 1986; Romani 
et al. 1986; Rickett & Coles 1988; see Appendix A of FM93). They found that the rms differential 
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refractive angular wander is 

(A9 2 ) 1 / 2 = 

Here r sep is the separation between pierce points through the scattering medium. For a dou- 
ble source (like B2 2050+36) with a component angular separation of 9 scp , r sep = D9 scp . The 
functions F(a) and G(a) are normalization functions that differ from group to group, but typi- 
cally by no more than a factor of a few. For reference, Cordes et al. (1986) predict (A9 2 ) 1 / 2 = 
0.18 mas (-Dkpc/A C m) -1 ^ 6 ^^ 3 for a Kolmogorov spectrum in a screen located Z?kpc kpc distant, pro- 
ducing an angular broadening of Oj mas, and observed at a wavelength of A cm cm. 

In general if a > 4 refractive effects like angular wander dominate over diffractive effects like 
angular broadening. Thus one would expect to see the apparent position of a source wander by 
more than the width of its scattering disk. If a < 4, the converse is true. 

Table 4 summarizes the 0.61 MHz measurements of the angular separation between the two 
components of B2 2050+36. We have estimated (A6* 2 ) 1 / 2 using two different methods. Hence- 
forth we shall also distinguish between the ensemble average expected differential angular wander, 
(AO 2 ) 1 / 2 , and the quantity estimated from the separation measurements, {Of.) 1 / 2 . First, we calcu- 
lated (0 2 ) 1 / 2 us ing the standard estimator for the variance and using the three values in Table 4. 
This estimate is (9 2 ) 1 / 2 = 2.3 mas. Second, we estimated (9 2 ) 1 / 2 to be comparable to the maxi- 
mum difference between component separations at the different epochs, max(A#). This estimate 
is max(A#) = 4 mas. 

Figure 3 shows (A9 2 ) 1 / 2 as a function of the density spectral index a, as predicted by Cordes 
et al. (1986), Romani et al. (1986), and Rickett & Coles (1988). We also plot our observational 
limits (9 2 ) 1 / 2 , 3(9 2 ) 1 / 2 , and max(A#). In producing Figure 3 we have assumed D = 2 kpc. We have 
considered (A9 2 ) 1 / 2 as a function of both D and a. However, the dependence on D is fairly weak, 
(A6* 2 ) 1 / 2 oc D-V6 for a = 11/3. Reasonable ranges for D produce little difference in (AO 2 ) 1 / 2 . 
The value D = 2 kpc is motivated by the distance to the Cygnus superbubble (Bochkarev & Sitnik 
1985), through which the line of sight to B2 2050+36 passes. 

Figure 3 shows clearly that along this line of sight the upper limits on (AO 2 ) 1 / 2 constrain 
a < 4. The largest allowed value, a = 3.86 occurs if we adopt an upper limit of 3(9 2 ) 1 / 2 and 
use the expression for (A9 2 ) 1 / 2 from Rickett h Coles (1988). Using a more stringent limit on the 
differential angular wander, max(A#), the expressions for (A9 2 ) 1 / 2 from Cordes et al. (1986) or 
Romani et al. (1986) we find a < 3.7. We adopt a = 11/3, as a limit based on (9 2 ) 1 / 2 we consider 
to be overly strict. 

The standard estimator for the variance, which we have used in determining (9 2 ) 1 / 2 , assumes 
that the measurements are independent. In this case the independence of the various measurements 
can be assessed by comparing the interval between the different measurement epochs with the 



/ \ (a-4)/2 

\G(a)(^y t 



a < 4; 
a > 4. 



(5) 
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Table 3. Best-Fit Diffractive Scattering Model Parameters for Eastern Component at 1 GHz 



Parameter 


Adopted 
Value 


Range 


0i 


3.9 mas 


3.5-4.3 mas 


Os 


7.0 mas 


6.98-7.03 mas 


a 


3.52 


3.51-3.53 


Derived Scattering Parameters 



SM a 10- 2 - 1 kpc m - 20 / 3 
h < 2000 km 



a Kolmogorov density spectrum, a = 11/3, as- 
sumed. 



Table 4. B2 2050+36 Component Separation Measurements 



Epoch 


$scp 


Reference 




(mas) 




1984 August 


60 ± 3 


Mutel & Hodges 1986 


1988 June 


56 ± 2 


FM93 


2000 July 


56.1 ± 1.0 


this work 


Differential Angular Wander Upper Limits 




2.3 mas 




max(A6>) 


4 mas 
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Fig. 3. — Differential angular wander (A^) 1 / 2 as a function of the density spectral index a. Three 
expressions for (A0 2 ) 1 / 2 are shown, that of Cordes et al. (1986) (dot-dash line), Romani et al. (1986) 
(dotted line), and Rickett & Coles (1988) (dot-dot-dot-dash line). All expressions are taken from 
Appendix A of FM93. The abrupt change in slope at a = 4 occurs because no attempt has been 
made to join smoothly the expressions for a < 4 and a > 4. The solid horizontal lines represent our 
upper limits on the differential angular wander of B2 2050+36 as given by the standard variance 
estimator of (0 2 ) 1 / 2 and 3(0 2 ) 1 / 2 . The dashed horizontal line represents our upper limit on the 
differential angular wander as given by max(A#). 
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refractive time scale, 

At r = ± = 20yr ( -?-) ( V , (6) 

where v is the bulk velocity of the material. FM93 assumed v = 50-100 km s -1 , for which At r = 2- 
4 yr. In this case, all three measurement epochs (Table 4) could be considered independent, though 
the first two only marginally so. However, if the characteristic velocity of the scattering medium 
is smaller, by as little as a factor of 2, the refractive time scale would be long enough that the first 
two epochs (1984 August and 1988 June) could not be considered independent. 

If the characteristic velocity is as small as 10 km s -1 , then less than one refractive time 
has elapsed since the first observation and none of these observations can be considered to be 
independent. Even if this was the case, we would expect to see some change over the 12 yr time 
span separating the second and third epochs (1988 June and 2000 July). That essentially no change 
in the angular separation is seen indicates that max(A#) provides a reasonably stringent upper limit 
on(A# r 2 > 1/2 . 

3.2.2. Refractive Intensity Fluctuations 

In concert with refractive wandering, the large-scale density fluctuations will also produce 
intensity fluctuations as more or fewer rays are "steered" into the line of sight. In the case of 
B2 2050+36, this will produce independent variations in the flux density of its two components. 

FM93 compared the flux densities of the two components between their observations and those 
of Mutel & Hodges (1986). They found that the flux density of the components was consistent 
with no change in the 4 yr between 1984 August and 1988 June. Table 5 summarizes these earlier 
measurements and those reported here. 

Because of changes in telescopes, receivers, correlators, and possible intrinsic changes, we 
do not compare the absolute flux densities of these earlier epochs to our measurements. Rather 
we compare the flux density ratio between the two components. As Table 5 shows, the western 



Table 5. Flux Density of B2 2050+36 at 0.61 GHz 



Epoch 


Eastern 


Western 


Ratio 


Reference 




(Jy) 


(Jy) 






1984 August 


3.8 


0.6 


0.16 


Mutel & Hodges 1986 


1988 June 


3.9 ± 0.1 


0.6 ± 0.1 


0.15 ± 0.03 


FM93 


2000 July 


2.864 ± 0.006 


0.374 ± 0.014 


0.13 ± 0.01 


this work 
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component has remained at a flux density of approximately 15% that of the eastern component for 
the past 15 yr. 

Goodman & Narayan (1985) present expressions for the rms refractive intensity fluctuation m r 
as a function of density spectral index a (see also Rickett 1986). In general, m r ~ 10% for a < 4 
and m r > 50% for a > 4. We determine an observed value m r ~ 1.5% from the flux density ratios 
in Table 5, suggesting a < 4. For a < 4, the strong scattering, i.e., large SM, determined from the 
observed scattering diameter (Table 3) means that the rms refractive intensity fluctuation will be 
only a slowly-changing function of a. From our observed value of m r we find (Goodman & Narayan 
1985, equation 3.1.5) a < 3.85, consistent with the value found by fitting the angular diameters 
(Table 3). 

3.2.3. Comparison with PSR B 2020+ 28 

The pulsar PSR B2020+28 (8?7 from B2 2050+36) has been the subject of many monitoring 
programs designed at least in part to measure refractive scattering along the line of sight to it. In 
all cases, it is found that diffractive scattering dominates refractive scattering, as we find for the 
line of sight to B2 2050+36. 

Graham Smith & Wright (1985) used dynamic spectra obtained at 0.41 MHz to measure 
the frequency drift rate and estimated that the diffractive scattering angle was at least 5 times 
that of the refractive scattering angle. Using a structure function analysis of a 0.43 GHz flux 
density monitoring program, LaBrecque, Rankin, & Cordes (1994) find diffractive scattering to 
dominate over refractive scattering. Bhat, Gupta, & Rao (1999) monitored the pulsar's flux density 
at 0.33 MHz and estimated the diffractive and refractive scattering angles toward this pulsar. They 
found that the diffractive scattering angle was roughly 10 times larger than the refractive scattering 
angle. They also estimated a density spectrum spectral index of a = 3.54. Stinebring & Condon 
(1990) were unable to observe any refractive modulations of the pulsar's flux density in their multi- 
frequency monitoring program, consistent with diffractive scattering dominating on this line of 
sight. 

4. Intrinsic Structure and Kinematics of B2 2050+36 

A key assumption in our analysis of the differential image wander is that there is little intrinsic 
change in the separation of the components. Given that on milliarcsecond scales, many sources 
show evidence of jets and bulk motion, this assumption clearly must be justified. 

We have compared our 1.67 GHz observations with those of Phillips & Mutel (1981). This 
comparison is somewhat difficult because, in general, different antennas at different locations were 
used in the two VLBI arrays. Thus, this comparison is necessarily crude, but it will be sufficient 
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to show that our assumption is reasonable. 

Phillips & Mutel (1981) cite a separation of 60 mas between the two components from obser- 
vations in 1979 September. If we fit our 1.67 GHz observations with only two components, we find 
a separation of 58 mas, less than that of Phillips & Mutel (1981). Even allowing for differences 
in arrays and uncertainties in the fitting, we conclude that any intrinsic change in the separation 
between the two components is probably less than 1 mas. Over the 20.8 yr separating these mea- 
surements, any intrinsic motion is probably less than 0.05 mas yr -1 . At a redshift of 0.354 for 
B2 2050+36 (de Vries et al. 2000) and for an assumed Hubble constant of 65 km s _1 Mpc" 1 , this 
proper motion corresponds to a projected linear separation velocity of c. Our limit is similar to 
that found by Shaffer, Kellermann, & Cornwell (1999) for CTD 93, a source with a similar compact 
double structure. 

5. Conclusions 

We have conducted a multi- frequency suite of observations on the source B2 2050+36. We 
have used these observations to confirm and extend previous conclusions, principally those of FM93, 
regarding interstellar scattering along this line of sight. 

Our 0.61 MHz observations form the third epoch of VLBI observations at this frequency. Our 
observations are separated by 12 to 16 years from the previous two epochs. As do previous obser- 
vations, we find the structure of B2 2050+36 at this frequency to be dominated by two components 
approximately 56 mas apart. Within the uncertainties of the various measurements, this separation 
has remained unchanged for the past 16 years. We can place a reasonably strict limit of 4 mas on 
any differential image wander caused by refractive interstellar scattering. In turn, this limit implies 
that the electron density power spectrum along this line of sight has a spectral index near the 
Kolmogorov value, with a value of 4 being highly unlikely. As the diffractive scattering diameter 
of the source at this frequency is 21 mas, we conclude that diffractive scattering dominates along 
this line of sight. 

The most significant source of uncertainty in this limit is probably the velocity at which the 
scattering material is drifting past the line of sight. A "high" velocity of 100 km s _1 would mean 
that the separations measured at the three epochs are independent. A "low" velocity of 10 km s _1 
would mean that none of the epochs are independent. Even if the velocity of the scattering material 
is low, the separation between our observations and the previous ones is probably long enough that, 
if refractive scattering were important on this line of sight, some change in the apparent separation 
of the components would be expected. Nonetheless, future observations of B2 2050+36 at 0.61 GHz 
or similar observations of other sources on lines of sight for which the refractive time scale could 
be estimated would be useful. 

Our observations are in good agreement with a variety of observations on the pulsar PSR B2020+28 
(line of sight 8?7 from B2 2050+36). These observations suggest that diffractive scattering also 
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dominates refractive scattering along the line of sight to the pulsar. 

Our quasi-simultaneous observations range from 0.33 to 8.4 GHz. We have used these multi- 
frequency measurements to constrain the frequency dependence of the angular diameter of one 
component of B2 2050+36. At frequencies below 1 GHz, the structure of the eastern component 
of B2 2050+36 is dominated by scattering. Fitting for both the density spectrum spectral index a 
and the scattering diameter Oj, we find a = 3.52 with the value of 4 again highly unlikely. At a 
nominal frequency of 1 GHz, we find the scattering diameter to be 0d = 7 mas. We also place a 
conservative upper limit on the density spectrum inner scale of l\ <C 2000 km. 

We also compared our 1.67 GHz observations to those of Phillips & Mutel (1981), acquired 
approximately 21 years previously. We find no evidence for any change in the separation at this 
frequency. This places a limit on the projected linear separation velocity of the two components 
of c. 

The National Radio Astronomy Observatory is a facility of the National Science Foundation 
operated under cooperative agreement by Associated Universities, Inc. Basic research in radio 
astronomy at the NRL is supported by the Office of Naval Research. 
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